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ABSTRACT 

Context. The presence of the lithium-6 isotope in some metal-poor stars is a matter of surprise because of the high values observed. 
Non-standard models of Big Bang nucleosynthesis and pre-Galactic cosmic ray fusion and spallation have been proposed to explain 
these values. However, the observations of this light isotope are challenging which may make some detections disputable. 
Aims. The goal was to determine ^Li/^Li for a sample of metal-poor stars; three of them have been previously studied and the 
remaining two are new for this type of study. The purpose was to increase, if possible, the number of lithium-6 detections and to 
confirm previously published results. 

Methods. Spectra of the resonance doublet line of neutral lithium Li i 670.78 nm were taken with the High Dispersion Spectrograph 
at the Subaru 8.2 m-telescope for a sample of five metal-poor stars (-3.12 < [Fe/H] < -2.19). The contribution of lithium-6 to the 
total observed line profile was estimated from the ID-LTE analysis of the line asymmetry. 

Results. Observed asymmetries could be reproduced assuming isotopic abundance ratios ^Li/'Li of the order of: 0.004 for BD 
-^26° 3578, ~ 0.010 for BD +02° 3375 and G 64-37, 0.025 for BD +20° 3603 and 0.047 for BD -04° 3208. We found that these 
results were very sensitive to several of the assumptions made in the analysis, in particular, the treatment of the residual structure in 
the analysed spectra. Our final estimates for the errors are respectively A^Li/'Li = ±0.028, 0.029, 0.039, 0.025 and 0.039. 
Conclusions. The ^Li/'Li ratios for the sample are comparable to or even lower than these error values, so that detections of lithium-6 
can not safely be claimed despite of the high resolving power (R ~ 95 000) and S /N (400-600). 

Key words. Line: profiles - Stars: abundances - Stars: atmospheres - Stars: Population II - Galaxy: evolution 
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> 1. Introduction 



The abundances of light elements in metal-poor stars are a 
good in dicator of what hap pened early in the Universe and the 
Galaxy dAlpher et al. I [1941) . The origins of the two stable iso- 
topes of the light element lithium, ^Li and ^Li, are still puz- 
zling. Their derived abundances in metal-poor stars are difficult 
to reconcile simultaneously with standard predictions. Different 
production sites in the early Ga laxy have been discussed: Big 
Bang nucleosyn thesis (Wagoner et alJI 19671: iNoUett et alj|1997t 
ICoc etal .''2004': 'Jeda mzikii2q08.) . cosmic ray spallation and/or 
a-fusion (Reeves etal. IVangioni-Flam et al. 2000) and 

neutrino spallation (see e.g. Vangi oni-Flam et alj 11996 ). The 
most widely discussed picture is that lithium-6 has a cosmic 



ray-induced origin (princ i pally g- fusion) (Vangioni-Fl am et al 



200C : lAoki et 



[gin (princ ip: 
al] l2004t 



lAsplund et alT [2006c iRolUnde et al 



20061 and re ferences therein) and lithium-7 a predominantly pri- 



mordial one ( ISpite&Spitelll982l) . 

The observed abundances of ''Li are used to dis- 
crimi nate between differen t production sites: pre-Galactic 
(e.g. ISuzuki & Inoug l2002t iRollinde et af] [20 06. 200 ^, the 
Galactic interstellar medium (e.g. iMeneguzzi et alj 19711; 



Fields &01iv3 119 99^. superbubbles (e.g. iParizot & Drurvl 
1999 ; Vangioni-F lam et al. 20 0Q|), fast supernova eject (e.g. 
Fields et al. 1996; Nakamu ra et a n 120061) . in-situ stellar flares 
(Tatischef f & Thibaud 2007i) . etc. The problem with this iso- 
tope is that detec tions for metal-po or stars are very challeng- 
ing and very rare. I Smith et al] (Il993h reported the first probable 
detection of ''Li in HD 84937 (''Li/^ Li = 0.05 + 0.02) which 
was later confirmed by other works dHobbs & Thorburnlll994l ; 
iSmifli et alJlT998l;ICavrel et alj|l999l) . A few years after this first 
detection, a detection for anot her star, BD -i-26° 3578 (%i/'^U = 
0.05 +0.03. ISmith et al.]|l998[) . was rep orted. Po ssible detections 
for two other stars were provided by ICayreTe t al. (1999) and 
iNissen et all (l2000b : BD+42°2667 (''Li/^Li = 0.10) andG271- 
162 (^Li/^Li - 0.02 + 0.01) respectively. There are results pub- 
lished for other metal-poor stars also, however the values are 
mostly upper limits. Our group determined ''Li/^Li < 0.018 for 
the metal-poor star HD 140283 ( Aoki e t al. 2004); this st ar may 
have par tly depleted its initial content dCavrel et alj|l999 ). The 
work of Asplund et al.l (l2006 l) has the highest number of pub- 
lished detections in metal-poor stars: nine. The detections for 
these stars appear to lie in a plateau that, if confirmed, would be 
difficult to explain within the standard framework of fight ele- 
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ment production. One would expect an evolution of the lithium- 
6 abundances with metallicity because of the contribution of 
Galactic cosmic ray spallation. Analyses of ^Li for very metal- 
poor stars ([Fe/H] < -2.5) are still very limited, so new samples 
of a few stars have a large impact. 

A possible solution to the apparent existence of a plateau 
and also to the high abundance of that plateau could lie in a 
pre-Galactic origin o f the cosmic ra ys (Suzu ki & Inoue 2002; 
iRoUinde et all l2006h (see, however IEvoU et a l. 2008). Non- 
standard models of Big Bang nucleosynthesis (BBN), generally 
involving particles of uncertain character, are also considered 
as a possible solution. The importance of establishing whether 
or not ^Li has been reliably detected in metal-poor stars is thus 
readily apparent: the previously published detections have been 
a significant motivator of investigations of the role of supersym- 
metry in modifying nucleosynthesis during the Big Bang. The 
detection (or otherwise) of ''Li could therefore potentially pro- 
vide crucial, rare observational constraints on the nature of cer- 
tain hypothesised supersymmetric particles and dark matter can- 
didates, some of which cannot be studied experimenta lly with 
curre nt or near-term experimental facilities (e.g. Baill v et all 
l2008h . 

Although detections for ^Li are more common than for 
^Li, the origin of th e heavier isotope is also quite intriguing. 
ISpite& Spite! (11 9821) observed essentially the same lithium con- 
tent, the so called Spite value, in a sample of metal-poor dwarf 
stars with effective temperatures 5800 < T^ff < 6250, despite 
their metallicities ranging from -2.2 < [Fe/H] < -1.1. A plau- 
sible interpretation for the plateau is that it is close to the pri- 
mordial value. ^Li is the only light isotope besides H, D, ""He 
and "^He that, according to standard BBN models, is predicted to 
be produced in significant amounts in the Big Bang. Predictions 
for the primordial lithium content based on standard Big Bang 
nucleosynthesis (SBBN) in combination with the baryon den- 
sity inferred from the Wilkinson Microwave Anisotropy Probe 
Mission (WMAP) give a value two or three times high er than the 
observed plateau value, A ( Lifl = 2.65 vs 2.1-2 3 dex ( Coc et al.l 



l2004tlAsplund et al.l2006trHosford et al.l2009tlAoki et al.l2009 

and references therein). 

Lithium is destroyed by nuclear reactions at high tempera- 
tures. In the presence of mixing, material from the stellar at- 
mospheres can reach depths of high temperatures where lithium 
can burn or can mix with material that has burnt it. According to 
some stellar models, little or no li thium-7 destruction is e xpected 
for stars in the Spite plateau (e.g. lDelivannis et"ani 19901 and ref- 
erences therein) so that these stars should retain their initial con- 
tent. However atomic diffusion is expected to act in warm main- 
sequence star s and alter the initial photospheric chemical com- 
position (e.g. iRichard et al.ll2002, 2005 ). Stellar models includ- 
ing diffusion and/or stellar rotation have been suggested to rec- 
oncile the observations of metal-poor stars wit h the SBBN pre - 
dictions, e.g Pinsonneault et al. ( 1999) and Rich ard et al.l(l2005l) . 
Observers have searched for signatures of mixing by measur- 
ing the lithium-7 abundance dispersion and by studying the be- 
haviour of the abundances with stellar temperatures and metal- 
licity. There is not consensus yet whether the Spite plateau 
is a plateau or not and whether lithium abundances are dis- 
persed orjio^_(e^. Bonifacio & Molaro 1997; Rvan et al. 1999, 
20001; lAsplund et al..,2006; Bonifacio et al..,.2007,; .Hosford et al.i 



2009tfAokietal.ll2009h~ Lithium-6 is however more fragile than 
lithium-7 (burning at temperatures T ~ 1.6 x lO^K versus 
T ~ 2.0 X 10^ K), so an observable content of the lighter iso- 



A(Li) = logA'u/iVH + 12 



tope sets constraints on the degree of depletion of the other 
one. Stellar models without diffusion or rotational mixing pre- 
dict part destruction of lithium-6 in the pre-main-sequence 
phase for most of the stars with published lithiu r n-6 detections 
(iProffitt & Michaudll98 9^'DeUvannis et a l.lll990l:PVsplund et al.l 
|2006|) . According t o lAs plund et al. (2006) neither the rotational 
mixing models of Pinsonneault et al.l fll99 9) nor the diffusion 
models (including turbulent diffusion) from Richard et al. (200^ 
2005) seem to give a solution to both lithium problems simul- 
taneously. A high depletion of lithium-6 would be difficult to 
reconcile with observations of the other stable isotope because 
a very high initial content of lithium-6 would generally imply 
a significant lithium-7 production by cosmic ray nucleosynthe- 
sis. Nevertheless, rece nt results for extremely metal-poor stars 
(iBonifacio et al.l l2007l) show a correlation between lithium-7 
abundance and metallicity which may have an interpretation in 
terms of Galactic evolution by cosmic ray spallation. 

Given the importance of the published lithium-6 detections 
for the understanding of Big Bang and cosmic ray nucleosyn- 
thesis, it is paramount to confirm the detections and, if possible, 
to add more. This is particularly important because of the very 
weak signal presented by ''Li: a small depression of the red wing 
of the ^Li doublet, at a level comparable to possible sources of 
systematic error. We present here ^Li/^Li estimates based on ob- 
servations taken with the Subaru 8.2m-telescope for five metal- 
poor stars, of which three are in common with previous studies. 
The observations are described in § |2] while the stellar param- 
eter determinations and the spectral synthesis are in § |3] In § |4] 
we present the stellar and instrumental broadening calculations. 
In §|5]we explain the procedure followed to get our best isotopic 
abundance ratio estimates, and we present the results and their 
sensitivity to different uncertainties, while in §|6]we discuss the 
results in terms of the observational quality and the sensitivity to 
the spectral analysis. We finish with conclusions in §|2] 

2. Observations 

The ''Li/^Li estimates presented here are based on the spectral 
analysis of the observed Li i resonance doublet line at 670.78 nm. 
Spectra were taken for eight metal-poor stars with the High 
Dispersion Spectrograph (HPS. Hoguchi et"an 120021) mounted 
on the Subaru 8.2m-telescope on the nights of May 17-19, 2005 
(UT). Results for five stars are presented here; the observations 
of three other stars were not completed in this run to achieve the 
data quality required for this study. Values for exposure times 
are listed in Table[Tl together with the V magnitude, date of the 
observations, exposure time, the central wavelength of the in- 
strumental setting employed and the measured heliocentric ra- 
dial velocity of the star. Total exposures were split into several 
shorter ones for cosmic ray rejection. Three of the stars were 
observed in two consecutive nights. Observations were taken us- 
ing a 2048 X 4100 LEV CCD ( 13.5yL(m pixel size), a sHt width 
of 0.4arcsec (0.20 mm) and two standard instrumental settings: 
one with central wavelength at 608. 1 1 nm (StdYd) and the other 
at 620.78 nm (StdYb) (0.002 nm/pixel dispersion). Special care 
was taken for the calibration; a high number of flat-field frames 
(~ 20) were acquired and ThAr lamp spectra were taken at sev- 
eral times. Fast rotator stars were also observed, in case their 
spectra were required for flat-fielding correction. 

The observational data were reduced under iraf using the 
echelle package and following the standard procedure: bias sub- 
traction, flat-fielding correction, spectra extraction, wavelength 
calibration, normalisation and combination of individual spectra 
(spectral orders were not merged). For the wavelength calibra- 
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tion, an average of 20-40 Th lines per spectral order were iden- 
tified. The wavelengths of these lines were used to convert pixel 
values into wavelengths. The wavelength solutions are based on 
a fifth order polynomial fit to the transformation. Typical values 
for the rms of these fits are 0.06 pm. If necessary, the fits were 
interpolated in time before being applied to the observed spec- 
tra. The normalisation consisted of dividing the spectra by their 
cubic spline fit with a wavelength scale of 0.3 nm. Spectra were 
co-added and weighted according to their exposure times. 

The co-added spectra are of very high quality. The signal- 
to-noise ratio (S IN) per pixel was estimated from the measure- 
ments of noise on both sides of Li i 670.78 nm line in a window 
of 0.07 nm, values are in the range 417-61 1. The resolving power 
{R - /I/A/Ifwhm) was estimated from the full width at half max- 
imum (FWHM) of a Gaussian fit to a set of ThAr lines lying 
in the spectral order that contains the lithium resonance doublet 
Une; results are in the range R -90000-100000. 



3. Stellar parameters, model atmospheres and 
spectral synthesis 

The observed spectra were analysed using spectral synthesis. 
The first step in the analysis was to get an estimate of the fun- 
damental parameters describing the stellar atmospheres: effec- 
tive temperature Teff, the surface gravity logg, the metallicity 
[Fe/H] and the m icroturbulence (t- The va lues listed in Table[T] 
were taken from Nissen et al.l (l2002l 12004) . with the exception 
of the metallicity. The choice of the stellar parameters is not 
very imp ortant for the accuracy of the ^Li/^Li final results (e.g. 
ISmith et al. 1993; Asplund et al. 2006), however we have tried to 
have the maximum consistency. Although for some stars there 
are more recent values than the values assumed here, we have 
chosen those which were determined using the same indicators 
and in a similar way: {b - y) and (V - K) colors for Tetf, and 
Hipparcos parallaxes for \ogg. The uncertainty in the measured 
par allax of G 64- 3 7 is v ery high so an estimate coming from 
the ISchuster et al.l (1200 4*) calibration of the Stromgren indices 
based on Hipparcos parallaxes was used instead. Our metallic- 
ity estimates come from the measured equivalent widths of a 
set of Pen lines: ^614.92, ^624.76, /1 641.69, ^643.27, and 
A 645.64 nm. Iron abundances were derived and taken as a metal- 
licity estimate. The solar values were determined from t he mea- 
sured equivalent widths of these lines as in Garcia Pere z et all 
(1200 6) [but now using new 3D solar C and O abundance values]. 
For G 64-37, the Fe ii lines above were not detected, but one at 
523.4 nm was, giving [Fe/H] - -2.97. As our metallicity calcu- 
lation for this star would be based on only one line, we preferred 
to take the value from Nissen et al. (2004) that is determined 
from several lines: [Fe/H] - -3.12. 

Equipped with the stellar parameters (see TablelTJ, model 
atmospheres were co mputed using the marcs code version in 
lAsplund et al.l ( Il997h . The code solves the equations of hy- 
drostatic equilibrium and energy equilibrium assuming that the 
atmospheres have either plane-parallel or spherical symmetry 
(ID); we assumed plane-parallel geometry in our computations. 
The radiation and the matter are assumed to be coupled, so the 
radiation can be described by the local atmospheric conditions; 
radiation is said to be in local thermodynamical equilibrium 
(LTE). The same assumptions were made for the spectral syn- 
thesis carried out with the bsyn code from the Uppsala package. 

The validity of these assumptions is a subject of discussion. 
Atmospheric convective motions can produce asymmetries in 
spectral line profiles. These asymmetries cannot be reproduced 



with the use of ID model atmospheres, so in principle our results 
for isotopic abundance ratios based on observations of the pro- 
file asymmetry may be in error. A more appropriate treatment 
of the convection, like the one employed in the construction of 
hydrodynamic model atmospheres (3D), may be required for a 
bett er estimate of the li ne asymmetry due to the isotope lithium- 
6 (Asplu nd et al.ll2066l) . However Asplund et al. investigated the 
effects of the use of 3D model atmospheres on their ^Li/^Li cal- 
culations and found no significant differences with respect to the 
use of ID models. The investigation was done assuming LTE 
conditions. Studies o f the formation of the Lii 670.7 line (e.g. 
ICarlsson et al.l [T994 ) suggest that lithium abundances derived 
from this line observed in spectra of metal-poor stars may be 
affecte d by departures from LTE by up to 0.03 dex. lCavrel et al.l 
(120071) showed that for the metal-poor star HD 74000, ID-LTE 
and 3D-NLTE results for ^Li/^Li disagree (0.024 vs 0.006). The 
3D-NLTE analysis of Cayrel et al. suggests there is no ''Li in the 
star. One of the purposes of our work is to compare our results 
with values reported in other works which are basically ID-LTE 
results; that is why we have made a ID-LTE assumption. The 
3D-NLTE effects cannot be neglected when discussing implica- 
tions of our results for the modelling of chemical evolution of 
the universe. 

The Li i 670.7 line is a multiple! and has two isotopic com- 
ponents, ''Li and ^Li. The assumed values of wavelengths and 
\oggf of these transitions are the same as in Smith et al. ( 199^. 
Laboratory measurements suggest they are known with a high 
accuracy. Atomic data for the Fe ii hnes used for metallicity de- 
terminations were tak en from lNissen et al.l (|2002|) and from vald 
(iPiskunov et al.|[T995l) for lines other than these. 

4. Stellar and instrumental broadening 

Our estimates for the isotopic abundance ratios (^Li/^Li) are 
based on the comparison of the observed and synthesised profile 
of the resonance Li 1 670.78 nm doublet line. Before the compar- 
ison, the synthetic spectrum has to be convolved by a function 
that reproduces the instrumental and the stellar line broadening; 
for this work we have assumed a Gaussian function. We adopted 
the approach of including the stellar and instrumental broaden- 
ing contribution in a single free parameter, the broadening pa- 
rameter Abroad (= A/ipwHivi [kms"']), that was determined from 
the analysis of up to 26 observed spectral lines: Nai (588.9 nm), 
Mg I (470.2, 552.8 nm), Ca i (558.8, 612.2, 616.2, 643.9 nm), Ti ii 
(444.3, 446.8, 450. 1, 453.4, 456.3, 457.1 nm), Fe i (441.5, 452.8, 
491.8, 492.0, 523.3, 526.9, 532.8, 537.1, 539.7, 561.5 nm) and 
Fell (458.3, 492.3, 501.8 nm). 

Different R values for the different lines could compromise 
the accuracy of our determinations, however R measurements 
from ThAr lamp lines in the range A — 550 to 700 nm did not 
show significant systematic variations with wavelength (AT? ~ 
1%). The spread of values at a given wavelength was higher, 
though still not so significant (~ +5%). 

The other two free parameters of the analysis were the rest 
wavelength of the observed spectra and the chemical abundance 
(log ex) of the line studied |3 The free parameters were varied 
until the best fit to the observations was achieved, using the 
technique. 

Values of the broadening parameter versus the measured 
equivalent widths of the spectral lines are shown in Fig. [T] for all 
the observed stars of our sample. Not all lines were measurable 
in all stars. The data in Fig.[T] suggest that lines stronger than 

- logex = logA'x/A'H + 12 



4 A. E. Garcia Perez et al.: ^hifhi estimates for metal-poor stars 

Table 1. Observational information 



Star 


V 


Date 


Exp 


Time 


Inst. Setting 


Vmd H 






[Fe/H] 


^, 




(mag) 


(UT) 


(s) 


(nm) 


(kms ') 


(K) 


(cgs) 


(dex) 


(kms ') 


BD-04° 3208 


9.99 


18-05-05 


4x 


1800 


620.78 (StdYb) 


58.7 


6338 


4.00 


-2.28 


1.5 






19-05-05 


2x 


1800 














(j o4-i 1 


1 i. 14 


1 O AC AC 


7 X 


1800 


0/0. /8 (Std 1 b) 


OO 1 

6/. 1 


1 o 

631s 


A 1 

4. lo 


Q 1 O 

-3.1/ 


1.5 






19-05-05 


7 X 


1800 














BD+02° 3375 


9.93 


19-05-05 


4x 


1800 


620.78 (StdYb) 


-397.1 


5855 


4.16 


-2.19 


1.5 


BD+20° 3603 


9.69 


18-05-05 


1 X 


1800 


620.78 (StdYb) 


-241.6 


6092 


4.04 


-2.22 


1.5 








2x 


2400 


















19-05-05 


1 X 


1800 














BD+26° 3578 


9.37 


17-05-05 


3 X 


1800 


608.11 (StdYd) 


-128.2 


6239 


3.87 


-2.33 


1.5 








1 X 


2400 















~ 7 pm require lower broadening values so mean values from 
the weaker lines were determined for Abroad and taken as our best 
estimate, see Table|2] The weakest lines observed have equiva- 
lent widths in the range 0.4 to 3.8 pm depending on the star For 
comparison, the observed Li lines have equivalent widths in the 
range 1.6 to 3.1 pm (see TableO, values are shown in Fig. [T] as 
dashed vertical lines; only BD -1-20° 3603 has a lithium line that 
is much weaker than the minimum equivalent width included in 
its broadening calculations (2.6 vs 3.8 pm). We assumed that for 
this star the broadening parameter does not change significantly 
with equivalent width in the weak-line regime. 

The mean values of the broadening parameter, listed in 
Table|2] were well determined for all the stars with typical stan- 
dard deviation values of 0.15-0.40kms The observed Cai 
612.2 line and its best fit is shown in Figure|2] where you can 
appreciate the low sensitivity of the fit to changes in the broad- 
ening parameter of the order of ±0.5 kms ', this value is higher 
than the maximum of the cr^ values listed in column 3 in Table|2l 



5. Isotopic abundance ratios 

5.1. Calculations 

The ^Li isotope is shifted by ~ 0.015 nm relative to ^Li in the 
neutral lithium resonance doublet line. This produces an asym- 
metry in the already asymmetric profile of the red wing of the 
line. Under the assumptions of ID model atmospheres, observed 
asymmetries may only be due to the multiplet configuration of 
the lithium resonance line and to the isotopic shifts (for 3D ef- 
fects, see discussion in §[3]l. 

The best fit to the observations of the lithium line was de- 
termined from ?L analysis (x^ - J^iOi - S i)^ /cr^^^; where 
Oi - normalised observation in pixel /, 5, = normalised synthe- 
sis and (Ts/N - noise based on the S /N) with three free parame- 
ters. A first, approximate, continuum normalisation was done in 
iRAF. The continuum level of the spectra was redetermined and 
fi xed prior to the analysis, not taken as a free parameter as 
in lAsplund et al.l f2006) . The continuum levels were determined 
by linear interpolation between two windows (see FigureO, one 
on each side of the lithium line. Different regions were selected 
for different stars in order to minimise the presence of structure. 
The presence of ripples that remain after this continuum normal- 
isation compromises the accuracy of the continuum level as well 
as the determination of the S IN. We return to this issue below. 

With the broadening parameter and continuum level fixed, 
the free parameters of the analysis were the lithium-7 abun- 
dance A(^Li), the isotopic lithium abundance ratio ''Li/^Li and 
the wavelength zero point of the observations. The three free pa- 



i BD-04° 3208 
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Ti « 
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1 BD+02° 3375 



OO 



I BD+20°3603 



1 BD+26°3578 



o O 



5 

Eqw [pm] 



10 



Fig. 1. Broadening parameter versus the measured equivalent 
width of the different spectral lines used in the calculations. 
Lines of different atoms are shown with different symbols: Fe 
(circles), Ti (stars), Ca (triangles), Mg (squares) and Na (dia- 
monds). 

Table 2. Broadening parameter values and their standard devia- 
tion. 



Star 


4'broad 






(kms"') 


(km s-') 


BD -04° 3208 


6.61 


0.18 


G 64-37 


6.24 


0.39 


BD -1-02° 3375 


5.20 


0.32 


BD +20° 3603 


5.84 


0.23 


BD +26° 3578 


6.49 


0.16 



rameters were allowed to vary until x^ was minimised. They 
were changed in steps of: 0.02 dex in lithium abundance, 0.01 
in isotopic abundance ratios and 0.0002 nm in wavelength shift, 
with typically allxl3xll points grid. Wavelength points in a 
range from 670.743 to 670.823 nm were included in the;t'^ cal- 
culations (36-37 wavelength points). Our best ^Li/^Li estimates 
for all the stars are listed in Table[3] (column 8) together with the 
S IN around the resonance neutral lithium line (column 3), the 
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Table 3. Equivalent widths and S /N measurements of the Li 1 670.78 nm together with the best results for the;^f^ three free parameters 
and the sensitivity of the analysis. 



Star W^i SjN AA Ai'Li) "Li/"Li A^Li/'Li 

(pm) (nm) (dex) 













ncor 


cor 


best-esti 


cont 




Abroad - ( 


BD -04° 3208 


2.37 


602 


0.01036 


2.27 


0.047 


0.043 
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Fig. 2. Spectral observations of the Ca 1 612.2 line for the five ob- 
served metal-poor stars (circles) and spectral synthesis for three 
values of the broadening parameter, the best estimate (solid line) 
and the best estimate + 0.5 km/s (dashed lines). 



wavelength shift (column 4) and the lithium-7 abundance (col- 
umn 5). FigureslH |5] |6] |2]and [8]show the observed metal-poor 
spectra and their best synthetic fits for ""Li/^Li = 0.00, 0.05, 0. 10 
together with their residuals (observations - synthesis in %). The 
goodness of the final fits can be appreciated from the l^x^ values 
(Kx^ - X^- the minimum value) presented in these figures. Note 
that the values were not divided by the degree of freedom so 
the values are non-reduced values. 

5.2. Contamination 

Possible "contamination" of the lithium line was investigated by 
comparing the spectra of the observed metal-poor stars with the 
spectra of hot fast rotator stars taken the same night. This com- 



parison is shown in Figure[3] where two plots are shown for stars 
observed in two different nights, each corresponding to a differ- 
ent date of observation. The associated hot fast rotator spectra of 
the four metal-poor stars G 64-37, BD +02° 3375, BD +20° 3603 
and BD +26° 3578 show tiny features of peak-to-valley range 
~ 0.003-0.006 (in normalised flux units) that suggest that their 
observed lithium line suffer from some contamination. At most 
one or two of such features could be identified with diffuse in- 
terstellar bands. We presume the majority is residual fringing or 
residual flat-field eiTors. Note that, although the observed rip- 
ples would not significantly affect common abundance determi- 
nations, they can have a great impact on the isotopic abundance 
ratios. One approach to quantify the effect of these ripples in the 
''Li/^Li results is to divide the metal -poor spectra by the spectra 
of the hot fast rotator stars. The division has two effects: one is 
to alter the flux in the Li line directly, and the other is to require a 
revision of the continuum fit, which again, affects the normalised 
flux in the line. The correction is not perfect; there are features in 
the observed metal-poor spectra that do not match the features in 
the fast-rotator spectra (e.g. BD -04° 3208) so we cannot be sure 
the lithium line is free of ripples at the level of accuracy that we 
require for our determinations even after this correction. Results 
for the isotopic abundance ratios coming from non-corrected and 
the corrected (divided) spectra are presented in Table[3] columns 
6 and 7 respectively. 

The observations of BD +02° 3375 were from only one night 
and, unfortunately, the S /N of the fast rotator spectra for that 
night was not very high, S /N ~ 250. To avoid degrading the 
quality of the analysed spectra, the division by a fast rotator 
spectrum should be done only when necessary and so for this 
star it was not done. Had the division been done for this star, 
the signal-to-noise of the analysed observations would have de- 
creased {S/N ~ 437 to 300) but the *'Li/^Li results would 
have not changed significantly, 0.007 vs -0.001 (see values in 
Table[3]). For this star we have adopted the ''Li/^Li value based 
on the non-corrected spectra (higher S/N) as our best estimate, 
while the ripple effect (''Li/^Li|cor-^Li/^Li|„cor) is included in 
the error bar 

The problem of the low S /N of the fast rotator spectrum was 
not such a big issue for the other stars observed that night, G 
64-37, BD -04° 3208 and BD +20° 3603, as they were observed 
also on a second night. Therefore, the S /N per night required to 
achieve the high target S /N was lower than for the stars which 
were observed on only one night. Still, the low S/N makes it 
difficult to judge by eye whether there is, or not, contamination. 

Observations of G 64-37 suggest no contamination for the 
first night, whereas it is hard to say for the second night. The 
^Li/^Li corrections for G 64-37 are very small, A^Li/^Li ~ 
0.007, so for this star as well, no correction was done. The com- 
parison of the spectra of BD -04° 3208 and its associated fast 
rotator does not show any sign of strong contamination. Hence, 
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our isotopic ratio results coming from non-corrected and cor- 
rected spectra are quite similar, 0.047 vs 0.043, our best estimate 
is based on the non-corrected spectra. 

Results for BD -1-20° 3603 were more affected by the division 
of spectra: *Li/^Li ~ 0.004 (non-corrected) vs 0.025 (corrected). 
We consider it more appropriate for our best estimate to adopt a 
value based on the corrected spectra, given that there are indica- 
tions of contamination. 

This was also the case for BD -1-26° 3578; ripples are ob- 
served in the hot spectrum at the wavelengths of the Lii 
670.8 nm line, therefore we applied the correction to get 
our ^Li/^Li best estimate, 0.004 (corrected) vs 0.074 (non- 
corrected). The match between the ripples in the metal-poor 
spectrum and in the associated fast rotator is quite impressive 
(see Figure[3j. The observed residual fringing in the fast rota- 
tor spectra recorded on May 17*, 2005 was greater than on the 
other nights when a different instrumental setting was used. BD 
H-26° 3578 was the only programme star observed on that night 
(see Table[T]). The observed spectra shown in Figures|7]and[8]are 
the divided spectra. 

5.3. Sensitivity 

The sensitivities of our results to the different assumptions made 
in the analysis have been studied. Concerning the normalisation, 
changes of 0.1% in the continuum level, where the typical line 
depth is 10%, has an effect of 0.02 on the ^Li/^Li results, so the 
normalisation of the spectra is an important issue. 

After the initial continuum fit (see §|2l), the spectrum of BD 
-1-26° 3578 showed a remaining tilt, so it was re-normalised in 
iRAF using a linear function similar to that employed in the 
analysis. After this was done, the match between the ripples in 
the spectrum of BD -1-26° 3578 and in the associated fast rotator 
spectrum was excellent, hence we have used a linear fit in our 
analysis. For stars with only minor differences in the measured 
continuum level on either side of the lithium line, the use of the 
mean value of the continuum fluxes, rather than a linear interpo- 
lation, does not have a big effect on the ^Li/^Li results. However, 
for BD +26° 3578 and BD +02° 23375, the isotope ratio changed 
by A^Li/^Li = 0.005 when using the mean-value approach. 

The level measurements can change with the spectral win- 
dow choice, different areas were selected for different stars and 
they are shown by horizontal bars in Fig. [3] A change of choice 
within the range of possibilities would affect the results signif- 
icantly for all the stars except BD +20° 3603. Even a shift of 
the window in the left wing by < 0.1 nm will change the iso- 
tope abundance ratios significantly, the changes ranging from 
A''Li/^Li = -0.037 (BD -04° 3208) to 0.028 (BD +02° 3375), 
see column 9 in Tabled BD -04° 3208 and BD +02° 3375 are 
the stars with the highest uncertainty in the continuum location 
among our sample. For these two stars, the use of the fast rotator 
spectrum did not help for discriminating between two different 
windows in the left wing, either because the small structures ob- 
served in the fast rotator spectra did not match the ones in the 
metal-poor spectra or because the spectra of the comparison star 
had a low S /N. The result of this was that two normalisation lev- 
els were possible. Summarising, isotopic results for these stars 
have a high uncertainty because of the uncertain continuum lo- 
cation. 

Another aspect to take into account is the number of wave- 
length points included in the calculation, in particular in the 
red wings of the lithium line. A decrease of typically four or five 
wavelength points (~ 0.010 nm) in the red wing affects the iso- 
tope ratio results by a few tenths of one percent for three stars 



(Tabled column 10). The stars G 64-37 and BD +26° 3578 were 
more affected by these changes, A''Li/^Li - 0.023 and 0.014 re- 
spectively. 

Results also show sensitivities to the assumed values of 
macroturbulence; we determined new isotope ratio results for 
broadening parameters values decreased by their cr^. We find a 
mean sensitivity of A^Li/^Li ~ 0.01 (see values in column 11 
of Table[3]l, although the sensitivity for BD +26° 3578 is even 
higher, 0.016. 

We have investigated the effects of changing the stellar pa- 
rameters for one particular case, BD -04° 3208, according to 
their uncertainty: Areff = 100 K, Mogg = 0.2, A[Fe/H] = 
0.2 dex and A^t = 0.2 km s"'. Isotopic ratio determinations are 
not very sensitive to these changes directly, A^Li/^Li < 0.004 
(see top row of Table|4|. A second effect of changing the stel- 
lar parameters is to alter the inferred broadening parameter (see 
second row of Tabled. For example, an increase in the micro- 
turbulence results in lower Abroad values; the broadening parame- 
ter compensates particularly for those parameters most strongly 
connected with the profile width: A^t and \ogg. We have shown 
already (Table[3] column 11) that the isotopic ratio is quite sen- 
sitive to the broadening parameter The total effect of increas- 
ing the microturbulence will be a ^Li/^Li increase of ~ 0.003 
(fourth row of Table|4]l, -0.003 from microturbulence directly, 
and +0.006 from the resulting change in Abroad (third row of 
Table|4|l. A total error associated with the uncertainty in the stel- 
lar parameters including the broadening parameter is estimated 
from the quadrature sum of the values in fourth row in Table|4] 
and is A^Li/^Li ~ 0.006. Note that this excludes the direct effect 
of the uncertainty in the broadening parameter, which is given in 
Table|3] column 1 1 . Summarising, the two first rows of Table|4] 
list the partial sensitivities to the stellar parameters, the third row 
the partial sensitivity of ^Li/^Li results associated to the change 
in the broadening parameter as a result of changing the stellar 
parameters, and the last row the total combined effect (stellar 
parameters + broadening parameter). The quadrature sum of all 
the uncertainties are listed in the last column. 

A summary of all the individual sensitivities are given 
as A^Li/^Li values in Table|3] (columns 9 to 11) and the fi- 
nal row of Table|4] The quadrature sum of these plus the 
^Li/^Lilfor-^Li/^Lilnror values were taken as our best estimates 
of the total ^Li/^Li errors (Table[3] column 8). We considered the 
division by a hot fast rotator spectrum to be essential, not merely 
optional, for BD +26° 3578, so the ^Li/^Li error bar for this star 
does not include the uncertainty associated with that correction, 
otherwise the error bar would be overestimated. 

Another concern is whether asymmetrie s are intrinsic to even 
single lines. Following lCavrel et al.l (l2007h . we have studied the 
asymmetry of the observed profiles of two Fe i lines at 639.36 nm 
and 649.50 nm in the spectrum of BD -04° 3208. For that, the 
profiles were re-scaled to the same line depth, they were aver- 
aged and after that they were convolved with a Gaussian function 
of width 2.39 km s ' to corre ct for the different t hermal broaden- 
ing of the Fe and Li lines (see lCavrel et al.l2007l) . The magnitude 
of the asymmetries shown in the determined profile depends on 
the assumed normalisation; the presence of fringing can intro- 
duce asymmetries. Continuum levels that are acceptable choices 
can produce an asymmetry whose additional depth is of the same 
order of magnitude as that produced by a displaced line of equiv- 
alent width 4% of the main one. Therefore, despite the quality 
of our observations, it is difficult to determine whether the asym- 
metry is real or just a product of the normalisation. 
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Table 4. Sensitivity of the ^Li/^Li and the broadening parameter results to the stellar parameters uncertainties for BD -04° 3208. 
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5.4. Summary 

Our best estimates of the ^Li/^Li ratios in these stars are set out 
in Table[3] column 8 and displayed in Fig.|9l The tabulated un- 
certainty is a quadrature sum of the contribution from uncertain- 
ties in: whether or not division by a hot fast rotator is performed 
(§ l5.2l and Table[3] columns 6 and 7), ~ 0.010, changes resulting 
from different choices of continuum windows (§ |5.3| and Table|3] 
column 9), ~ 0.022; changes resulting from different choices 
of the wavelength range over which the observed and synthetic 
spectra are compared via the statistics (Table|3] column 10), 
~ 0.010; the star's broadening parameter (Table|2] column 3 and 
Table[3] column 1 1), ~ 0.012; and the star's atmospheric param- 
eters (TablelH row 1), ~ 0.006. We note, but do not include, fur- 
ther possible eiTor sources: whether a mean value of the contin- 
uum windows or a linear interpolation between continuum win- 
dows is adopted as "the continuum"; and whether even single 
spectral lines are asymmetric. The tabulated uncertainties range 
from +0.025 to +0.039, and are dominated by the uncertainty in 
the continuum, even though the S /N exceeds 400 in every case. 
The reason that the continuum placement continues to exert such 
an influence on the outcome of the analysis lies substantially in 
the problem of small structures seen in the reduced data, and 
that although they are small by most subjective measures, they 
are significant relative to the very weak signal of ^Li, the spec- 
tral line of which has a central depth of probably only < 0.003 
normalised flux units (i.e. < 4% of the depth of the ^Li, which 
in turn is ~ 0.08 flux units). Based on the extreme sensitivity 
of the isotope ratio results to alternative choices in the spectral 
analysis, we are certain that we could not infer ^Li/^Li to better 
than +0.025, or perhaps more realistically +0.040, despite hav- 
ing data of 7? = 95 000 and S /N ~ 400 to 600. Moreover, we 
note that the formal eiTors which we quote are not the maximum 
possible values and in some cases will underestimate the real 
error, particulary since they result from the quadrature sum of 
several distinct contributions. The continuum from cr^ in partic- 
ular is only a Icr value. 

6. Discussion 

According to models of big bang nucleosynthesis and cosmic- 
ray spallation, low lithium-6 abundances are predicted in stars 
of low metallicity. Detections of this isotope in these stars would 
have thus a strong impact on production models of light ele- 
ments. All the tests we have done, e.g. comparing the results 
based on non-corrected and corrected spectra, on diff'erent nor- 
malisation levels, on different wavelength ranges covered in the 
analysis, suggest that our best estimates for ^Li/''Li are uncertain 
to at least 0.025-0.040. Given than our best estimates of ''Li/^Li 
are of the same order of magnitude as these values, we do not 
claim any significant detection. 

One of the goals of this paper was to observe targets whose 
lithium-6 content was already determined from observations 
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Fig, 4, Top: the observed spectrum of the Li 1 670.78 nm line for 
BD -04° 3208 (circles) along with synthetic spectra for three 
diff^erent ''Li/^Li values: 0.00 (soHd), 0.05 (dotted) and 0.10 
(dashed). Middle: differences between observations and spectral 
synthesis given in %. Bottom: non-reduced (three free param- 
eters, 36 wavelength points) variations with ''Li/^Li; the values 
are given relative to the minimum (Ax^). 



taken with other spectrographs, and show whether, or not, the 
published results were corroborated by our observations. There 
are three objects in common with other studies: BD +26° 3578, 
BD -04° 3208 and BD +20° 3603 (color symbols in Figure|9l 
see the electronic edition of the Journal for a color version of 
the figure). 
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Fig. 5. Same as in Figure|4]but for G 64-37 and 37 wavelength 
points. 

The first was observed in ISmith et aP (Il998h and 
lAsplund et alj (l2006l) (red symbols in Fig.|9]l. They re- 
ported different isotopic ratios results: *Li/^Li - 0.05 + 0.03 
(S/N ~ 420) and ^Li /^Li = 0.01 + 0.013 {S /N ~ 680) respec- 
tively. lAsplund et al.l s uggested th at the differences between 
their value and that of ISmith et aU could be due to the differ- 
ences in the determination of the broadening parameter. We 
have used a similar method to that employed by A splund et al.L 
the difference being that we did not use a rotational profile 
(metal-poor stars have low r otational velocities). We obtain 
Abroad = 6.50 vs 5.65 km S-' ( lAsplund et al.1 l2006l ). We have 
data of similar S /N {~ 642) but slightly lower resolving power 
{R ~ 95 000 vs 1 15 000). Our best estimate for this star (^LifU 
= 0.0 04 + 0.028) agrees with the result from Asplund et al.l 
( 120061) . Much higher values ^Li/^Li = 0.074 would be deter- 
mined from the non-corrected spectra, but based on the match 
between the residual features seen in our spectra of this star 
and the corresponding fast rotator (see Figure|3] bottom right), 
we are confident that division (correction) by the fast rotator is 
essential for BD -1-26° 3578. Our observed bright fast rotator star 
spectra suggest a residual fringing of the order of 0.4 - 0.6% 
(peak-to-valley) whereas Asplund et al. report a value much 
lower for theirs < 0.2% (+0.1%). The amplitude of the fringe 
patterns of the HDS spectra (peak-to-valley difference) at 



Fig. 6. Same as in Figure|5]but for BD -h02° 3375. 

670.0 nm is typically 5%, which results in a pattern of a few 
tenth percent after flat-fielding. The width of the fringe pattern 
is much larger than the Li absorption feature, and the effect 
on the line analysis is usually not severe. We note, however, 
that the effect is not negligible in t he setting ad opted on May 
17th 2005 (see §5.2). According to ISmith et al.L their residual 
fringing is negligible although they do not quote any magnitude 
so it is difficult to judge whether fringing can account for the 
differences between the results. 

Our b est estimate for BD - 04° 3208, the other star in com- 
mon with lAsplund et al.ll2066l (blue symbols in Fig.|9] see the 
electronic edition of the Journal for a color version of the figure) 
is 0.047+0.039 which agrees wefl with their value 0.042+0.019. 
Our error bar is higher so we do not claim any detection in con- 
trast to them. Concerning the quality of the spectra, we have 
slightly higher S /N (~ 600 vs 530). 

The other star in common with ' Smith et al.1 (1 19981) is BD 
+20° 3603 (green symbols in Fig.|9]l. We have data of similar 
quality to theirs (S /N ~ 430), though their R was a bit higher 
(~ 1 15 000). Our best estimate of 0.025 + 0.025 is higher than 
their best estimate of 0.00 + 0.02 but consistent within the uncer- 
tainties. If no fast rotation correction was done to our spectrum, 
a ^Li/^Li value closer to 0.00 would better reproduce the obser- 
vations. Differences in the assumed stellar parameters for these 
stars are not significant so they should not affect the results. 
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Stellar parameters assumed for the stars BD -04° 3208 and 
BD +26° 3578 are similar to one another The stars also have 
spectra of similiar quality. Their isotope abundance ratio results 
are different but agree within their error bars, 0.047 ± 0.039 vs. 
0.004 ± 0.028; the continuum location for BD -04° 3208 in par- 
ticular is quite uncertain. A change in the continuum level within 
the uncertaintity for this star would bring the best ^Li/^Li esti- 
mates for the two stars closer, ~ 0.01 vs ~ 0.00, though given 
that their temperatures differ by ~ 100 K, they may have experi - 
enced different degrees of *Li depletion (Delivannis et al.iri990h . 
The stellar parameters values assumed for BD +20° 3603 and 
BD +26° 3578 are slightly different and as in the previous case, 
their isotopic estimates are different but agree within their error 
bars. The agreement gets better if no correction to the observed 
spectra is made for BD +20° 3603. According to our error esti- 
mates, the ^Li/^Li ratios of these two stars are the most accurate. 

We are not aware of any other lithium-6 determinations in 
the literature for the two targets BD +02° 3375 and G 64-37. 
Unfortunately, the data quality for these two stars were among 
the lowest of the observed sample. Our best estimate suggests 
that we can only determine an upper limit of ~ 0.03 (Icr) for 
BD +02° 3375. The temperature assumed for this star is the 
lowest amon g our sample and low enough, according to the 
isochrones in lDehvannis et all ([l990), for the star to have suf- 
fered from severe ^Li depletion (~ 95% depletion). In fact, pub- 



lished detections of ^Li at metallicities similar to that of BD 
+02° 3375 are for hotter stars. Concerning G 64-37, the most 
metal-poor star in our sample, w e could estimate onl y an up- 
per limit of ''Li/^Li ~ 0.04 (Icr). lAsplund et al.l l2006' reported 
a ''Li detection for other star of [Fe/H] < -2.5, LP 815-43 
(''Li/^Li= 0.046 + 0.022). 

The ^Li/''Li values do not correlate with S /N; we deter- 
mined high isotopic ratio values from both high S /N spectra (BD 
-04° 3208) and low (BD +20° 3603), and the same applies 
to low values of the isotopic ratio. 

The principle conclusion of this study concerns not the iso- 
tope ratio itself but rather the uncertainties associated with the 
measurements. Despite having data of high resolving power 
(R = 95 000) and S/N ~ 400 to 600, a thorough examination 
of error contributions has led us to conclude that, for our dataset, 
we could not confidently determine the isotope fraction to bet- 
ter than 2.5% to 4.0%. Moreover, the detections achieved by 
other groups for metal-poor stars are around this limit, as we 
show with the ^Li contribution functions in Figure|9] determined 
from the linear fit to the ^Li observed abundances. It is quite 
possible that other studies have achieved a higher precision than 
was achieved in our study, despite the similarity of the initial 
datasets. Although, notice that Asplund et al. 2006 observations 
show lower levels of fringing than ours. In contrast, we interpret 
our results as indicating that extreme caution is required in fi- 
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Fig. 9. Derived logarithmic abundances of ^Li and *Li versus 
stellar metallicities. The linear fit to ^Li abundances is presented 
(solid line) along with ^Li contributions (dashed l ines). Our data 
are sho wn (stars) together with the LTE data from lAsplund et af] 
(I2OO6I) (sq uares), for the stars in common wit h our study also the 
data from lSmith etalJ (1 1998 *) (triangles) and 'Nisse n etal.lll999l 
(circles) are displayed. The filled symbols indicate detections, 
the bars are 3 cr upper limits and the non-black colors denote 
stars in common. [See the electronic edition of the Journal for a 
color version of the figure.] 



nally quantifying the uncertainties, in order to establish ^Li/^Li 
ratios to the level of precision required to be confident of detec- 
tions of ^Li at the 4% level or lower. We are not confident that 
we have significant detections of *Li in any of our five stars. 

7. Conclusions 

The observed spectra of the Li i resonance line were analysed 
for five metal-poor stars. Our best fits to the observations sug- 
gest that only the star BD -04° 3208 may have a high lithium-6 
content, ''Li/^Li = 0.047 + 0.039. However, the uncertainties of 
the calculations are comparable to the calculations themselves, 
such that positive detections cannot be claimed for any of them. 
In general the analysis is very sensitive to the assumptions made 
for the continuum location, the residual fringing treatment and 
even the wavelength range covered in the analysis. Our error bars 
(A'^Li/^Li ~ 0.025-0.040) tell us that detections may have been 
possible if the stars had a higher lithium-6 content (^Li/^Li > 
0.050-0.080) following a 2cr criteria. 

Not surpris ing, we obse r ved a symmetries in the profile of 
the Fe i lines as ICavrel et al.l (|2007b did. However, we have ma- 
jor concerns because the magnitude of the asymmetries depends 
on the spectral normalisation carried out, and the presence of 
residual fringing could also introduce spurious asymmetries. 

Our ^lAplA results are comparable to recent values in the lit- 
erature, but our error estimates are higher, the result of a detailed 
analysis of a range of uncertain aspects of the analysis associated 
with the significant presence of structure in our reduced data-set 



(summarised in § 15.41 ). In contrast. [Asplund et al.ll2006l tabulated 
only the uncertainties due to finite S /N — which they find usu- 
ally to be their dominant error — but not the wide range of error 
sources encountered in our work. This puts limitations on the 
detections so that only upper limits of ^Li/^Li were inferred for 
the stars observed in this work. 

One possible, but unjustified, generalisation of our results 
would be to question the reliability of all *Li/^Li measurements 
at levels of 4% or 5%. Making this generalisation would then 
call into question the evidence for a possible upper envelope or 
plateau in the ''Li abundance, which has caught the interest of 
astrophysicists in recent years. We believe there is a value in 
continuing research into the ^Li/^Li ratios in metal -poor stars, 
including the fullest treatment of possible error sources, to de- 
termine whether such an upper envelope or plateau exists. 
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